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Abstract 

We discuss the probability of random knotting for a model of self-avoiding polygons whose segments are 
given by cylinders of unit length with radius r. We show numerically that the characteristic length of 
random knotting is roughly approximated by an exponential function of the chain thickness r. 
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1 Introduction 

Knotted ring polymers such as knotted DNA 
molecules are synthesized in various experiments in 
chemistry and biology Jl], |[ || . The question on the 
topological constraint of the macromolecules was 
first formulated by Delbriick, Frisch and Wasser- 
man in the '60s. H || Then, topological properties 
of ring polymers have been studied through numer- 
ical simulations. §, 0, |, |, 0, (0], || || © @ 
Furthermore, the fractions of knotted species syn- 
thesized after the random cyclization of the circular 
DNAs are experimentally measured, and they give 
an experimental estimate on the effective diameter 
of the DNA molecules surrounded by the counter 
ions. [jl7| |l8); it is estimated so that the experi- 
mental values of fractions of some knots are consis- 
tent with the theoretical values of the fractions ob- 
tained for the hedgehog model |ll| of self-avoiding 
polygons. The model gives self-avoiding polygons 
consisting of cylindrical segments, where the length 
and the radius of the cylindrical segments corre- 
spond to the Kuhn length and the effective radius, 
respectively. 

Let us introduce the probability of random knot- 
ting. First, we assume an algorithm of constructing 
random polygons (or self-avoiding polygons) of N 
polygonal nodes. Then, for a given knot AT, we de- 
fine knotting probability Pk (N) by the probability 



of observing a random polygon (or a self-avoiding 
polygon) of the knot K. When K is trivial (K = 0), 
we denote it by Po{N). It has been found [pi |l3| 
that the probability of unknot has a decreasing- 
exponential dependence on the nodes A^: 



P (N) = C exp(-N/N ). 



(1) 



For the different models of random polygons and 
self-avoiding polygons, knotting probabilities of 
some nontrivial knots have been evaluated numeri- 
cally. || [I]], llll Through the simulations with 
the Vassiliev-type knot invariants, it has been found 
that for a given knot type K, the probability Pk(N) 
as a function of the number A~ of nodes can be well 
described by the following formula |lj, 



Pk{N) = C K 



exp (-N/N 



(2) 



From the numerical results, it is found that the pa- 
rameter Nk should have almost the same value for 
any knot. [[l4|, [l^] The result is also consistent with 
the simulation of self-avoiding polygons on the cu- 
bic lattice. |l(| Thus, we may assume that Nk 
should be given by a constant N c for any knot K: 
Nk = N c . We call N c the characteristic length of 
random knotting. 

In this paper, we discuss how the characteristic 
length N c depends on the excluded-volume param- 
eter of self-avoiding polygons. We introduce a new 
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method for constructing self-avoiding polygons with 
cylindrical segments with radius r. The method is 
a variant of the dimerization algorithm. It seems 
that the algorithm is more efficient than that of the 
hedgehog model. Then, we study the characteris- 
tic length N c (r) as a function of the cylinder radius 
r for the self-avoiding polygons. We give the esti- 
mates of the characteristic length by applying the 
formula (|l|) to the numerical values of the probabil- 
ity of unknot. Here the fitting parameters Co(r), 
No(r) are determined by the least square method. 

The numerical study on the characteristic length 
of random knotting as a function of the thickness 
parameter of polymer chains should be important 
in the study of knotted DNAs where the effective 
diameter can be changed according to the concen- 
tration of the counter ions. Furthermore, it could 
be also important in other systems of polymer rings. 
The length and the radius of the cylindrical seg- 
ments of the model correspond to the two fun- 
damental parameters of polymer chains, i.e., the 
Kuhn length (the statistical length) and the stiff- 
ness (or the excluded- volume parameter), respec- 
tively. In association with the numerical result 
of the rod-bead model |l3), Nechaev and Grosberg 
have addressed a conjecture |^(| that the charac- 
teristic length No as a function of some thickness 
parameter d should be given by the following: 



N (d) = N o {0)exp(d/t). 



(3) 



Here I denotes the Kuhn length. As we shall see 
later, however, the conjecture (H|) does not hold for 
the rod-bead model. On the other hand, the numer- 
ical simulation in this paper shows that the expo- 
nential behavior (||) gives a rough approximation for 
the new cylinder model of self-avoiding polygons. 

The paper is organized as follows. In §2 we intro- 
duce the cylinder model and explain the method of 
our numerical experiment of random knotting. In 
§3 we show the numerical data of our simulations 
and discuss possible relations between the charac- 
teristic length and the cylinder radius. We also dis- 
cuss for the rod-bead model how the characteristic 
length depends on the bead radius. Finally, we give 
some discussion in 54. 



2 The methods of simulations 

2.1 New method for constructing 
self-avoiding polygons of cylin- 
drical segments 

Let us explain the new method for constructing self- 
avoiding polygons of cylindrical segments. In this 
method, each segment of a polygon consists of a 
cylinder of unit length with radius r. The main 
structure of the algorithm is given by the follow- 
ing: (1) we generate a set of chains with cylindrical 



segments by the dimerization method; (2) we con- 
struct polygons by connecting two cylindrical self- 
avoiding chains with the concatenating method of 
Ref. |lCt| . In this method, we also calculate the 
statistical weight related to the probability of suc- 
cessful concatenation. 

Let us define the condition of an overlap be- 
tween a given pair of segments for our model of self- 
avoiding polygons. First, we assume that there is 
no overlap between any pair of adjacent segments. 
Second, for a pair of two segments which are not 
next-neighboring to each other, we assume that the 
two segments have no overlap if and only if the dis- 
tance between the central axes of the two cylinders 
is larger than 2r. Here we have considered that a 
central axis of a cylinder is given by the line seg- 
ment between the centers of the upper and lower 
disks of the cylinder. 

In the step (1) of the algorithm, we construct a 
chain with cylindrical segments randomly, and then 
we check whether there is an overlap or not for all 
unadjacent pairs of cylinders of the chain. If there 
is an overlap, then we give up the chain and con- 
struct a new chain from the beginning. If there is 
no overlap, then we keep the chain in the computer 
memory. 

The new model of self-avoiding polygons pro- 
duces free-joint rings with thickness parameter r, 
systematically. We can construct very long self- 
avoiding polygons with the cylinder radius r. It is 
known that the thickness of polymer chains plays an 
important role in the study of stiff polymers such as 
DNAs. [^l], ^2], ^3) We may assume that negatively 
charged DNA molecules with surrounding counter 
ions can be approximated by impermeable cylinders 
with the effective radius given by the screening ef- 
fect. Thus, the new algorithm can be applied to the 
study of knotted closed DNAs. 

2.2 The method for simulation of 
random knotting 

Let us describe the algorithm of our numerical sim- 
ulation of random knotting. For a given number 
N of polygonal nodes, we construct M polygons. 
Here the number M should be very large, such as 
M = 10 4 . Then, for a fixed knot type K, we enu- 
merate the number Mk of such polygons out of the 
M polygons that have the same set of values of the 
knot invariants with the knot K. We estimate the 
probability Pr- (TV) of knot K by the following 



Pk{N) 



M k 
M 



(4) 



We employ two knot invariants, the determinant 
Ak{— 1) of knot and the Vassiliev-type invariant 
V2 (K) of the second degree, as the tool for detecting 
the knot type of a given polygon, pll E5[ The values 
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Table 1: Values of the determinant of knot | 
&k(— 1) I and the second Vassiliev invariant i>g(K) 
for some simple knots. 
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Figure 1 : Probability Pq (N) of unknot versus num- 
ber N of nodes for the cylinder model: the nu- 
merical estimates of Pq(N) for r = 0.01, 0.03 and 
0.05 are shown by black circles, black triangles, and 
black diamonds, respectively, with error bars given 
by their standard deviations. Number N of nodes 
are given by lOj + 1 with j = 2, 3, . . . , 15. 
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Table 2: Numerical estimates of the characteristic 
length No(r) for the cylindrical model 
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In Fig. 1, the estimates of the probability of un- 
knot for the cylinder model with r = 0.01, 0.03, 
0.05 are shown against the number N of polygo- 
nal nodes. The error bars are given by the stan- 
dard deviations. For the cylinder model, we eval- 
uate the variance of the values of the probability 
Po(N) of unknot by taking into account both the 
statistical fluctuation of the number and that 
of the statistical weight |l(| appearing in the ring- 
dimerization procedure. 

From Fig. 1, we may confirm the exponential 
decay of Pq(N) for the cylinder model. We note 
that the result is consistent with all the numerical 
simulations for other models such as the Gaussian 
dynamical model Q , the Gaussian random polygon 
0. |l5) , the _hedgehog model [f[T| , and the rod-bead 



model j 13, |15J. The lines in Fig. 1 are theoretical 
curves given by the formula (|l]). The estimates of 
the parameters Co and N (r) for the ten different 
values of the radius r from r = to r = 0.9 are 



given m Table § 



of the invariants for some typical knots are given in 
Table @. 

3 Characteristic lengths of 
random knotting 

3.1 The exponential decay of the 
probability of unknot 

Let us discuss the numerical results of our simu- 
lations. Employing the cylinder model introduced 
in §2, we construct self-avoiding polygons with TV 
polygonal nodes where the cylinder radius is given 
by r. For a given number iV and a given value of 
r, we construct four sets of 10 4 polygons. [M — 

4 x 10 4 .) We consider for the polygonal nodes N 
fourteen different numbers from 20 to 150, and for 
the radius r 10 different values from 0.0 to 0.09. 



3.2 The characteristic length of the 
cylinder model 

A semilogarithmic plot of the numerical estimates 
of the characteristic length N (r) is given in Fig. 
2 against the cylinder radius r for the ten different 
values from r = 0.0 to r — 0.9. From Fig. 2, we see 
that the characteristic length Nq (r) as a function of 
the cylinder radius r can be roughly approximated 
by an increasing exponential function of r. Let us 
consider a fitting formula in the following 

N (r) = iVo(0) exp (ar) . (5) 

Here, there are two parameters to fit, i.e., N o (0) 
and a. Applying the formula (|^) to the numerical 
estimates of the characteristic lengths, we have the 
fitting line drawn in Fig. 2. The estimates of the 
fitting parameters are given by Nq(0) — 292 ±5 and 
a = 43.5 ± 0.6. From the viewpoint of the \ 2 test, 
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Figure 2: Characteristic length AT (r) versus cylin- 
der radius r for the cylinder model. The numerical 
estimates of Nq (r) listed in Table 2 are depicted by 
black circles together with their errors. The fitting 
line of eq. (5) is determined by the least square 
method with x 2 =42. 



Figure 3: Fitting curve of eq. (6) to the data of 
No(r) versus radius r for the cylinder model. The 
three fitting parameters of eq. (6) are determined 
by the least square method, which gives \ 2 = 2-2. 
The data points and their error bars are given in 
Table 2. 
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however, the curve does not fit well to the data. We 
note that \ 2 = 42 while the number of data points 
is ten with the two fitting parameters. Let us now 
consider another fitting formula 



N (r) = N (0)exp{f3r u ) 



(6) 



where there are three parameters No(Q), (3 and v. 
The new fitting curve is shown in Fig. 3 to the 
same numerical data points with Fig. 2. The best 
estimates of the fitting parameters are given by 
7V (0) = 271 ± 6, (3 = 29 ± 2 and v = 0.85 ± 0.02. It 
seems that the curve fits well to the data, although 
the value of \ 2 might be a little too small: \ 2 — 2-2. 
Thus, we may conclude that for the cylinder model, 
the characteristic length of the probability of un- 
knot as a function of the cylinder radius is given by 
the formula (|gJ) . 

We have evaluated the characteristic length 
No(r) for r < 0.1. For stiff chains such as DNAs, 
the radius r can be rather small such as r = 0.003. 
jy| We thus consider that the range: 0.0 < r < 0.1 
covers most important cases. Furthermore, when 
r > 0.1, numerical evaluation of the characteris- 
tic length becomes technically more difficult. We 
recall that even for the case of r = 0.09, the esti- 
mates No(r) is of the order of 10 4 , which is much 
larger than the polygonal nodes (N < 160) of our 
numerical simulations. 



3.3 The characteristic length of the 
rod-bead model 

Let us now discuss the characteristic length for the 
rod-bead model. Here we assume that the param- 
eter r denotes the bead radius. (Fig A) In Fig. 4, 
for the rod-bead model, the numerical estimates of 
the characteristic length Nq are plotted against the 
bead radius r. For an illustration, the estimates of 
the characteristic length No for the cylinder model 
are also shown in Fig. 4 against the cylinder ra- 
dius r. We note that when r = 0, the characteris- 
tic lengths of the two models should coincide. The 
two algorithms are essentially equivalent for r = 0. 



Table 3: Numerical estimates of the characteristic 
length of unknot for the rod-bead model. The val- 
ues for r = 0.0, 0.01, 0.03, 0.05, are calculated in 
this paper, while those of r = 0.1, 0.15, 0.2 are 
given in Ref. [15], which are also consistent with 
Ref. [13]. All the data are obtained by construct- 
ing 10 4 polygons (M = 10 4 ). 
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Figure 4: Two series of estimates of the charac- 
teristic length No(r) are shown by black triangles 
for the rod-bead model and by black circles for the 
cylinder model, respectively. The estimates for the 
cylinder model is given in Table 2, and those of the 
rod-bead model in Table 3. The two estimates for 
the two models coincide for r = 0. 
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From Fig. 4, we see that the numerical estimates 
of the two models really coincide for r — 0, which 
gives a partial consistency check of our numerical 
simulations. 

Nechaev and Grosberg discussed a conjecture 
that the characteristic length No as a function of 
the "chain thickness" d should be given by an ex- 
ponential form: N (d) = N (0) exp(27d/£), where I 
is the Kuhn length. |20j It seems that their conjec- 
ture is based on the numerical result |l| of the rod- 
bead model. However, we see from Fig. 4 that the 
characteristic length No(r) for the rod-bead model 
is not given by an exponential function of the bead 
radius r with 0.0 < r < 0.2. Let us consider a part 
of the range such as 0.1 < r < 0.2, and apply for 
the rod-bead model an exponential approximation: 
N Q {r) = 7Vo(0)exp(7r). Tneil; we find that thc 

parameter 7 is roughly given by 15. Thus, for the 
rod-bead model, the conjecture of the exponential 
dependence could be valid only for a limited range 
of the bead radius. Furthermore, it seems that the 
estimate of thc coefficient 7 is not consistent with 
the conjectured value of 27. 

4 Discussion 

Finally, wc discuss a possible connection of the 
cylinder model of this paper to the hedgehog model. 
Klenin et.al. |ll| studied the probability of unknot 
for the hedgehog model of ring polymers. The al- 
gorithm of the model is given by the following: (1) 
generating a set of vectors of unit length with a 



common origin (a "hedgehog"); and (2) setting thc 
minimal distance between any two unadjacent seg- 
ments larger than the segment diameter. 

It seems that the hedgehog model and the cylin- 
der model introduced in this paper produce almost 
the same results for the self-avoiding polygons with 
clylinder radius r. It is our conjecture that the two 
models should give the same results, although the 
two algorithms themselves are quite different. For 
the probability of unknot, the two values obtained 
by the two different methods are rather close. How- 
ever, it is not easy to see whether they are exactly 
the same or not. We shall discuss much more precise 
comparison between the two models in our future 
publications. 
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Figure captions 

Fig. 1: Probability Po(N) of unknot versus number 
N of nodes for the cylinder model: the numerical 
estimates of P a (N) for r = 0.01, 0.03 and 0.05 are 
shown by black circles, black triangles, and black di- 
amonds, respectively, with error bars given by their 
standard deviations. Number N of nodes are given 
by lOj + 1 with j = 2,3,..., 15. 

Fig. 2: Characteristic length N (r) versus cylinder 
radius r for the cylinder model. The numerical es- 
timates of No(r) listed in Table 2 are depicted by 
black circles together with their errors. The fitting 
line of eq. (5) is determined by the least square 
method with x 2 =42. 

Fig. 3: Fitting curve of eq. (6) to the data of N (r) 
versus radius r for the cylinder model. The three 
fitting parameters of eq. (6) are determined by the 
least square method, which gives \ 2 = 2.2. The 
data points and their error bars are given in Table 
2. 

Fig. 4: Two series of estimates of the characteristic 
length N (r) are shown by black triangles for the 
rod-bead model and by black circles for the cylinder 
model, respectively. The estimates for the cylinder 
model is given in Table 2, and those of the rod- 
bead model in Table 3. The two estimates for the 
two models coincide for r = 0. 
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